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ABSTRACT. 


The determination. af elastic: properties: oficrystals :by 
sound velocity measurements: is: facilitated by use-of 7a: 
perturbation solution to the: equations of motion: for:sound 
waves propagating In. the crystals... Transformation: equations 
can relate the salution for arbitrary porpagation. direction 
to an established crystal coordinate. system. A-method ‘of 
Gauguletinge the acoustic components: of the energy flow vector 
is developed for several high symmetry crystal_classes. The 
walidity of this method. is herein. confirmed and:-results are 
presented which are consistent: with experimental_data-.. A 
modern desk calculator was: employed: to determine- solutions 


for hexagonal, cubic and tetragonal’ crystals.. 
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A major portion of the theoretical development: presented. 
faeconis paper Wes previcusly been presented. by, Froi.) Jonnie 
Neighbours, Physics Department,. U.. &S.. Naval Postgraduate: 
School. The most recent paper, co-authored by Associate: 
Professor G. E. Schacher,. Physics Department,. U..S2.Navai_ 
Postgraduate School, is titled "Determination. off Elastic: 
Constants from Sound Velacity. Measurements: in. Crystals: of: 
General Symmetry." In addition,. several. investigators: have 
published methods similar ta this development... These:papers 


are included in the List of References... 





LL. SReCUer Oy 


The task of relating velocity measurements to elastic: 
constants of a crystal is complicated by the fact: that. the 
Beesenvectors which describe the direction CP yibrea on eed 
response to the sound wave are neither perpendicular nor 
parallel to the propagation direction. The energy flow may,. 
Mmewever, be described by modes of vibration which are quasi-- 
longitudinal (essentially parallel to the propagation. direc-- 
tion). In practice one vibrational mode dominates the: others: 
j/meroa Cliven propagation direction. In addition the transducer - 
used to measure the velocity has greater response to one- modé: 
of vibration so it is possible to separate modes by appropriate: 
Bmoice of transducers. 

Mer certain high symmetry erystal classes specmal direc— 
Pees Of propagation exist in which the respomse is pure 
longitudinal and pure transverse. Velocity measurements: along: 
these directions can be directly related to crystal elastic. 
@emsGants. However the energy flow is not necessarily paradiel. 
to the propagation. 

For either high or low symmetry cases the elastic: constants: 
@en be determined by solution of the crystal vibration equation 
of motion. This process may be simplified by appropriate 
choice of coordinate system and application of several assump-- 
tions. The resultant secular equation tis solved by a 
perturbation method and the values for the elastic: constants: 


may be refined by iteration. The constants thus obtained are 





linear combinations af the elastic: constants: defined -_by- 
conventional coordinate: systems: and: can. be=converted iby: 
Metacional transiormataoy.. Ihe azgile- of deviag1on on amie 
energy flow vector from the propagation direetion. is:of: 
moe cular. Tnterest im order, Gorconstruct azpnysical picguure 
of the crystal response.. This: angle: cam be: obtained:by 
application of the general energy fllow equations: to -:the.: 


coordinate system defined by the- propagation. direction. . 





II. SECULAR EQUATION 
The general form of Hooke's Law is 


Bogen OS aval (1) 


Sa 
Two common conventions of subscript notation for the above 


equation are shown below. 





FIGURE 1 


The strain, S, is related to displacement by the equation 


eee 
Suy Fe OU, FY) 
where u is displacment. In contracted notation Hooke's Law 
may be written as 


Wo SG. 8 (3) 





ie t= 
Y= 4 = Tag = Te = Ty (4) 
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The existence of a strain energy function. reduces: the: 36. 
eiiestic constants to 21 independent. constants: for the=most. 
penmeral case. This number cam be further: reduced: for: specific 
erystal calsses because of symmetry considerations... 

If an arbitrary coordinate system is: chosen. such that. the 
x' axis lies along the direction. af propagation the: equations 
motion cam be written in terms of the elastic~ constants. in 
the arbitrary coordinate system and displacement... In this: 
G@evelopment u’, v", w" are displacements: along the: x',.y*,. 
and z' axes respectively. The primed <eeeiton eee eee 
fact that choice of coordinate system ts: arbitrary:. The- 


equations of motion are 
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If plane wave propagation ts: assumed: the: partiai_deriva- 
tives with respect ta y" and z" vanish: and: the: equation. 
becomes 

ae Lae oon 
F = cr cone) tor, 2 V Psi & Wee. (6) 


x! apes eee ais: 








samtjlariy, the equations of motion. for the: remaining-directions 








are 
a... ess De 
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a w! 53 Orme Q- — om Wee 
Pos = @ ms = CH = +0? -+-C: ae 
2 44° re 5% | DD ye 


me agdilvronal assumption. ise that, the: solutions: are periodic 


of the form 


ew" = A. cos (kx'-at) 


J8 


oe A, cos: (kx'-wt) (7) 


wr = Ag cos (kx'-wt) 


The assumed periodic solutions are inserted in the equations 
of motion and the partial derivatives: are evaluated. The 
resultant equations may be further simplified by the relation 
V = w/k where V represents velocity and is not to be confused 
with v', the displacement in the y' direction. The equations 


may now be written 
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mae secular equation has three real roots corresponding to 
one guasi-longitudinal and two transverse modes: of: vibration. 
If the solutions are assumed ta he essentially: longitudinal 
@r transverse with small deviations from these: directions 


the secular equation. yields @ perturbation solution: : 
Th eee _ AWASC1E 
c. = \ rEg 


al ain ae Gea. 
— Cie) An(Chy -- Ch) 





eve that the amplitudes ofthe assumed periodic: selutions do 
mee aftect the calculation... The remaining soiutions:to the 


perturbation equations are 





| 
es | 
Cee = oV5 ce 26 Voss we Wahiatvote) 6-26 3 
Cy = © | 
comme OG 
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Cc! Cc? 
_ | _ 15 56 
Che pV2 7 fo + LL (9) 


In particular, if the terms in square brackets vanish the 
solutions correspond to pure mode response and. the elastic 
constants are determined directly from the measured velocities. 
When the square bracketed terms’ are present velocity measure- 
ments along several directions yield an adequate number of 


equations to determine all six elastic eonstants of the wsecii: 


10 





equation to first order. hese values are (he iis ee eee 
the square bracketed terms and the: solution is:refined. The 
primed elastic constants are linear combinations:of:those 
defined by the conventional coordinate: system.. The-rotational 
transformation relating the two systems: is: deveioped-in 
meetion IV. 

Whe amplitudes of the- pertodic: solutions: may be:determined 
from the equations of motioy.. If the amplitude:of:one mode 
fermcercst is fixed at an arbitrary vadue- the: remaining two 
amplitudes can be expressed in. terms: of. the: elastic.constants 
meee Teasured velocities... Kor example,.if: the- amplitude 


Ay £5 set equal to unity: the longitudinal: modé-solution is 


a. — CL. Ce 
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The components of the 
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The x’ component of this v. 


ENERGY FLOW 
Energy Flow vector: are?- 
+ w'X? 
+ wryt (iy) 


pee ails A 
Zr. 


ector can be calculated: as:follows: 
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Note that the plane wave assumption has been applied to yield 


only those terms which contain. partial derivatives;:with respect 


to x'. The transverse components of the Energy Flow vector are 
t a t avi ' OW! 
Ey ut ¢ Cle soy + Che ger t Che SET) 
er au! dv! ee chien 
a VOT ima "26 3 yg Wes sce (13) 
- au! re owl 
+ owt (Gly ayy t+ Cle ger t Che ger) 
bn BE. H. Love, A Treatise on the Mathematical Theory of 


Elasticity (New York, 19 


Gi. [s 
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aa oes SOS 
Eo u'r ¢ Cy 15 3x! us GG ox a Cl. 5 ie as 
dul av ow!’ 
+ A Ciy bear if SG ox i Nase Teer ) 
a ai ee) ON) tn 
shilpa \ Se Oxt 36 3x' tae ax! ) 
@m substitution of the periodic solutions. the x! component. 
becomes 
= ~  ¢ = Ga _)l Sie 
EL. wA,sin n fe! tae kA, sin n) # Cy 64 kA,sin. 1) 
Oe , Mpeihe. ~ . = 
+ Cis KA,Sin n)] + wA,sin. n [C4-.{-kA,. sin’ n) 
Cee (-kA, Sin nyet ¢! war kA sin Tel 
+ wA3sin n [ eee ] 
where n = kx! —- at. 
The above form may be further simplified by defining 
Ee 
en ma 
s Saree n 
The final form of the Energy Flow components is shown below. 
Pi, = Ci AS + CLAS + CLAS + 2C1-A,A, + 2C1,A,A, + 2C! 
xe Bd al aioe 2 Sis iS eoemt 115 yer eas 36" a S) 
= Cae + C!-AS + C!_ae 4+ Goin + Cy > AJA + (Cie + Ci AVA 
i yOrr 26 2 yo 3 66 i 2. 56 
65,9 Ge oe 


ILS 


Js 





2 t 2 t t t t 
5 + C2 5A2 + (Ce 6 + C3 AVA, 4 (Co, “+ C13)A,A2 


(Ch +: Cr6)ADAZ lS) 


The resultant Energy Flow vector determines deviation from 
Girection of propagation as shown in the diagram below. Note 


meae the angle 6 is the compliment of the polarization angle. 





FIGURE 2 
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IV. TRANSFORMATION OF ELASTIC CONSTANTS 


The desired elastic constants are dei imedmby yen tea 
unprimed coordinate system which is related to the crystal 
basis vectors a, b, ¢c in accordance with convention adopted 
Pe the 1.R.E. standards committee“ Specific examples of this 
convention are provided in Sections V, VI and VII. The primed 
coordinate system is chosen so that the x' axis lies in the 
meveplane. This leads to a relationship between the two 


Systems in terms of the direction cosines 


8 = sin 6 cos ¢ 
m= sin 6 sin ¢ (256) 
n= cos @ 





HIGURE 3 


“standard Committee, "Standards on Piezoelectric Crystals," 
Proc. I.R.E. 37 (September, 1949)> p. 1376. 


Ibs, 





The transformation. elementsi are given in: tapyiae eon 


below. 


ae, 





The elements have been reduced by use of: the: definition 


a? = ia a5 m°. 


For example: 


x = 2x! — me yt ~ena eZ! ' (13D) 


The transformation of elastic constants: may be:obtained 


from Hooke's Law;; 


Ts =e Ss. i 
ati ay a ( 9) 
If the transformattons off stress and strain: are: T! = BT, 
and S' = ¥S5 the elastic constants are related-by the 
expression 
‘—] 
' = 
where 4 = g and g ts the matrix shown-below. 
Pia m* né 2mn 2on 2em 
— 2*a° 0 @) 0 Lopmome 
eae merieos ac —2mn -22n emn“a° 
emma * Syria "= 0 QR a 338 Bees 
{a -m°ng ~ Na (12m )a— 2(1=2n°) 2s eal ie 
=, ln Lace 0 aie ee (2°=m- yo 


The Energy Flow equations and the rotational transformation 


may now be applied to specific crystal classes... 
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Ve. CRYSTALS OF HEXAGONAL SYMMETRY. 


In application ta certain. high symmetry crystai ciasses-: 
the method is simplified due to symmetry considerations -and-: 
the fact that certain primed elastic: constants: may vanish. 
mers teads to @ partially diagonal secular: equation: and-the: 
existence of semi-pure modes af vibration... For. hexagonal. 
crystals the I.R.E. convention. defines: the: unprimed: coordinate : 
Swevem such that the z-axis is parallel to: the: crystal_basis: 
meovwor ¢c, i.e. the sixfold symmetry axis;. and: the: x=axis:is: 
Pemeitel to a, the twa fold ratation axis... The: y-axis:is: 
Pemecendicular to zx and z to preserve: a right: handed: coordinate - 


system. With this definition. the: primed elastic: constants. 


become 
segilt ae Ce. ih. 
rae? pe! een 
salle — 
Gt. = sect. + na(I-2n°) (Co... + 20),) tna: 
15 il plese dy. 553 


ee “se = 9 


since C16 and Oey are zero the existence of a semi-pure mode 
1S implied. Physically this means that: the. vibrational response 


of the crystal lies along one of the primed system axes :for: 


ay 





one mode (pure transverse) and the remaining modes are quasi- 
longitudinal and quasi-transverse. The relative orientation 
@r the elgenvector components to the propagation ts shown 


below. 







¥3 - ( QF) 


Vo - Pure transverse 


FIGURE 4 


Note that if Cis were to vanish the secular equation 
would become diagonal and the response would be pure longi- 
@eerrnal or transverse. This situation may occur for materials 
pessessing hexagonal symmetry. .The existence of this acciden- 
tal pure mode was predicted by Borgnis2 for the case of 
Propagation in a direction that makes an angle 8 with the c 
axis defined by 

tan® 6 = oa a ae (23) 

11 13 Wy 

Mime above requirement corresponds to the fact that Oo Vanusmes 
and the secular equation is diagonal in the perturbation 
eelucion. 


Equations (22) are sufficient to calculate the velocities 


of each of the modes from known values of conventional elastic 


3n, E. Borgnis, "Specific Directions of Longitudinal Wave 
aan in Anisotropic icetae Phys. Rev. 98 (August, 
955), p. LOO. 
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constants. The elastic constants of numerous: hexagonal _ 
crystals have already been. determined ant: the: application of: 
the perturbation method may be demonstrated by comparison: of: 
calculated mode velocities with measured values:.. The: solilu-- 
meron for energy flow, however, requires: nine additional. 

7 ' ih T 1 fe as “1 1 | Bees 
constants i.e. Cio ee Coy os C35»: C5G>: Cane: 936s: Shes: Sug: 
woeme OL these may varnish as a result: of symmetry... For 

'" rT ¥! = rts : es ates : , 
hexagonal symmetry Crys C56» C. go and ©. he: are adi identically 


feo, The rémagining primed constants are 


, Be a 
Ci5 = @ C15 +n Ci3 
a 4. oes h Giese 
137 (1 -— 2n-a Cr + no (Ch 5 + Ca. NCin) 
C55 = na(C,. - Cj5) Cau) 
1 = he Bi 
C35 = -na(n Cay + (I-2n")Cy.~ — 0 C3 + 2(1- =n cmp 
Of. = SS nate. - co. = 2r,,,) 
6 2 11 ee KH. 


In addition, the fact that Clq and Cy. are- zero: leads:to- 
Ssamplified forms of the amplitude equations.. With" these 
eeuipaatications the equations for energy flow in hexagonal - 


crystals are 


: 2 a ee 
Pyr = CLAD + CheAS + CtaAg + 201 .A,A, 
Pye = (Cig + Ca)AAy + (Che + CheAoA, (25) 
Pp. = @Ylam + celae + (CE EO yA, A, + Cj AS- 
Zi Bis He 3> 3 Si, aS ae es 46! des 
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Some of the terms in the above equations: will vanish: when 
certain roots of the secular equation. ares applied: since-the 
amplitudes of the assumed. pertodic response: vanish:. In. the 
notation ee. below the second subscript: denotes’ the:root’ of 
mire secular equation that is used in the calculation.. For 


instance, A is the amplitude when. the root. V,. is- used. 


fail 
The amplitudes are 


| 8 pW gE O85. 
A= 1 yp =O Ay = 
iS 
An, = 0 Boy = Lo Aggy = 0: (26) 
oVy a 
Ay = _ a a 


The energy flow relationships: for the modes: of? vibrational 
response in hexagonal materials are 
For the root Vy (Quasi-Longitudinal) 


Zs el 


= " : 4: tT , | 
Pen Srpfar * Sgehtag F Maye haaisy- 
al 
BP. = CtA>. + ctlon@ a= (CN 42 C1. )A__A- 
Zin lise Jue bis 4 oul Sr) oe} hr Sale 


For the root V, (Pure Transverse) 


_ 2 

Pet CégAz2 

Ps = 0 (28) 
ME ne. 

Por = Chehss 
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For the roat Wes 


a az 
= " : 4 ’ = 
Pot ~ rity © Cheha + Shae yeh ss: 


Pee Go (29) 


ae: a Pe of | 
Se ee) ae ae ee oe ee 


(Quasi-Transverse:) 


rd 
| 


For this symmetry the above three equations: (27,.28;.and 29) 
Sew that the y" component of the energy flow vector.is 
always zero. Thus the angle yp defined on. page 14-is-:always 
zero and any deviation of the: energy filow vector. from the 
propagation vector takes place: in. the: x'-z2' plane:. 

in order to demonstrate the: validity. of? the=perturbation 
method solutions were computed: on. a Hewlett. Packard:9810A 
Geek @ealculator and compared to published: data: for: hexagonal 
SGaeowels,. specifically zinc amd beryl,. as: presented: by M. J. 
r. ene Transformation equations: using known..values of 
elastic constants were programmed and the: velocities obtained 
were in agreement to .05%.. The energy filow components were 
calculated and the angles of deviation fronr propagation 
Girecvion were consistent to within 2°° of those presented by 
Musgrave. The graphs on the following pages are examples of 
results obtained for x-z plane: propagation in hexagonal 
Crystals. In these graphs the deviation angle, 6,.is plotted 
feec lunctiom of the angie, ¢,. which the propagation vector 


makes with the z-axis for each of the three modes. In addition 


ure J. P.. Musgrave, "On the Propagation of Elastic Waves 


in Aeolotropic Media," Proc. Roy. Soc. 226 (1954) p..356. 


eile 





the sign of all deviation. angles: im. the: cubic: and: hexagonal 
eases has been reversed in. ordex to edmpare- Che resis eae 
Musgrave's data. Table I lists the known elastic: constants 


and densities required for the calculation. 


TABLE I 






Elastic Constants 
Ie 


Density 
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The metals Zinc and Magnesium are- two: different: examples 
of hexagonal material. The difference between.C,,-and C.., 
in the two materials affects the response: characteristics. 
In both materials the pues Tem eee mocellTs has: a:non-zero 
Geviation angie for energy flow... This: means: that. although 
the mode is a pure vibration, associated. with .the-root’ V, and 
iveme along the v' axis, the energy flow is. non=zero:and not 
coincident with the x' axis (the propagation direction). 

The response for Zinc indicates that no accidental.mode 
exists since the Quasi-Longitudinal (QL) response does not 
become zero except at the end points @ =0° (the ec direction) 


and @=90° (the x-y base plane). The absence of accidental 


OW. P. Mason, ed.,, Physdeal Acousties v.. Ills (New swore 
ieee Dp. 3c- 


Ce 





response can be verified by equation (23). The: Qtasi-- 
Transverse (QT) mode has @ zero value bet weeerthesenaeenntee 
but this does not imply pure mode response since: the: amplitude 
Ay 3 has a non-zero value at Ghis point.” PFor=Zinc’ (he Zeger 
effect of elastic anisotropy an the energy flow is: to: channel 
the energy of most of the waves toward the base:plane. 

An accidental pure mode ts observable in Magnesium. The 
Quasi-Longitudinal mode has a@ zero value at 6. ==51.6°. 
Analysis of the amplitude shows that. both Aj 3. and: A2,.are 
zero along this direction of propagation. Thus: the-amplitudes, 
and hence the vibrational responses. are collinear: with the | 
primed coordinate system. The Quasi-Transverse: mode:has 
Gesoentially the same form of response as Zinc: but. is:opposite-— 
ly directed. This is due to the difference in:.the-value of 


(C ~ C33) For tme twommateriters.. The magnitude: ol the 


i1 
deviation angles in Magnesium indicates that the-material is 


essentially isotropic. 


es 





(94d) 


f ' ’ : ) 
ONIZ NI 


HydOud 
Bee (ab mae Abs Hod yal de) [li ee. ae! | 
FIONVY NOILVOVdOd °“SA FISNY NOLZLYVIAAG 





°@ “SIXV-Z WOWMA NOILOAYId NOILVOVdOUd HO ATIONV 








Prrat lee 


NOLLV 





i 
rt 
: 


LON 


LY 


T 


LN 


C 


“YOLOAA NOILVOVdOUd WOWA NOTLVIAGA JO AIONV 


mi 


(ad) 


24 








(9qq) ‘9 *‘SIXV-Z WOUA NOILOGYIGC NOLTLVOVdOUd HO HTONY 
f ‘ : i 


} 


' NAISANDVW NI NoILvDVdoud 
oy SO, nA eee. _, Yoda, pt eel f i 
/ qIONY NOLTLVOVdOud “SA AIONV NOILVIANA 


1 





LN 


OL 


“YOLOUA NOLLVOVdOUd WOUMA NOTLVIANG. dO. dTONV 


es 


(ad) 


rau) 





Vi... CUBLC SYMMETRY Ch sii 


Calculations for cubic: symmetry. crystais: follow the same 
pecterm as that for hexagonal crysta@is:. The-conventienaz 
pees are pardiiel to the crystad basis: vectors:a; .b. and ec, 
@ee @lastic constants are. C,.5 = ee C74. 7 Cn3;3 and: Cee = Crys 
transformation af these: constants: leads: to: the-equations 


shown below. 


Cra = (gt +e + on. YC Be NS em © +-m@n©) (C15 + 2C ny) 
—— Co 2 An bc : . 

Cre = arma | age + 257 Cia ae + Cyy 

a e a —2 


Cle = na (2) +m —n “a ) (C}4- — Cjo- —- Cy) 
— =e 
ee 0 eee 1 


ees - 
Cfo = amma “(2° — mo) (0). — Co. — 2Cyy) 


For propagation in the x=+y,. x-z,.y=z:or. (110) planes the 
Censtants Clé and C56 vanish yielding a. partially diagonal 
Secular equation as for the hexagonal.case. The-velocity 
eameulations are tfdentical to. those for. hexagonal materiais. 


mee relations from the secular equations: are 













2 5 
6° ©6~—CCO GG /o 
we .- 2 (ct. +c) (Cry + or)? 4c! se (6)0 5 
fee ap Sag T Coe) = Ne ee 1ic55°° as 
C3ne) 


Further response analysis in. the. above planes yields 


amplitudes that vanish for certain modes of propagation and 
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elastic constants that vanish. by choice: of: plane:= of: propaga- 
moron. Thus jonly Che»: Clg 


the resultant energy flaw equations which are: similar:to. 


ane Co remain. to: be: inserted:in 


Eq. (25) for hexagonal symmetry.. The- graphs: on. the: following 
pages show solutions for Aluminum with sound: velocity’ propa- 
gation in the (100) plane (x-y plane) and the: (110) plane 

(6 = 45°). The results agree with Musgrave's:data.. The- 


elastic constants for Aluminum are 


C = 9.5 x rott dynes/em~ 
JIL 

Cis = Wei Se 10am dynes/cm~ 

Chy = 2.9 x, tott aynes/em™ 


The value for density is 2.65 gm/em-.. These: values:are:taken 
from Musgrave. 

The results for (100) plane propagation.in™= Aluminum: show 
the existence of the accidental mode at e = 45° which is~ 
presvicted by Borgnis.. ie (USE eptieers ot tne neseen eee 
curves is expected in the case of cubic: crystais:. Inraddition 
the deviation angle for the Pure Transverse mode-is-: zero. 
Thus the energy flow is colinear with the propagation direc- 
mena. For propagation in the (110) and (100) planes the 
angele, vw, 1s zero. 

The response for (110) propagation merits. particular 
attention. The Quasi-~Longitudinal mode deviation angle goes 
GO zero at 6 = 54.6° which is-the [111] direction. This is 
Mot an accidental mode but Is indicative, of: the face cma 


Cis is zero because of propagation direction. Thus the 
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secular equation. is dtagonal and. the three=modes:are-pure. 

In addition the velocities V5. and V. are. equal_and-the devi- 
MS 

ation angles for these roots: are: equal. and non-zero. The 


Boemomencn. of internal comical refraction exists:for tats 


propagation direction. It. cam be shown that PL, .and- 


2. yive 
ifas 


(De ton are constant.. Since the response:for A. is 
equal and opposite ta that. for’ A., the response:for: any 

Pepeat™ combination af the two-,. (making an angle: y,.with the 
y"-axis) would occur at twice: the angle y,.in: the:-y direction. 
In other words the energy flow describes a:cone:about the 
Prepeagation direction... It the: amplitude: for: one-mode is 
rotated the response: rotates’ in the opposite-.direction at 


twice the angular frequency of rotation. 
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VIL. ‘TETRAGONAL SYMMETRY. 


solutions for tetragonal crystals represent: a: step: 
forward since there is little data availabie=- to: be- compared 
mech. The L.R.E. conventions: define: the z-axis: parailet to. 
the basis vector c and the x-axis parallel to: the: basis: 
vector a, or the twofold symmetry axis.. The: two: types:of: 
tetragonal class materials are distinguishable: by the-fact- 
that Cié = 0 for Tetragonal I materials.. The: transformation. 
Or C€lastic constants is otherwise identical... The= transforma-- 


erenms for Tetragonal Il materials are 


Cii = (x" + m')Cy, a 22m Cy + 2hhaQ (Oza. + 20jy) tn 3. 
+ Ag nce. + Yam(2 — m)Cy¢ 
CEs = nea ete + m')C1, 1 ze mn a (C5 + Bee) 
- n°a°(20,, - Czz) + (1 - 2n®) Poy + 4- maa" (2°---m* Cx ¢ 
Che = 22°ma “(C,, — Cyy) HNC, Ha (e> — me) Cee. 
- Yama *(2° - m°)Cy¢ (32) 
Cis = ee ee + m')Cy4 - 29 °m°n oron na (1 -- 2n*)(C,. +-2C) 1) 
+ nal, _ 4emna "Ce - 4amna (2° —m°)C4¢ 
te = -ama (2° - m°)(C14 - Cyo — 2g) 
+ ree ne oan Cae 
Che = gmna-( 2° = m°)(Cy4 - Ci» AAG a2) 
- he “Ee Ai = 62°m~)C1¢ 


Sal 





Once again the choice of (100) orm (110) plane: propagation 
eliminates the constants ore and C1-¢ and. the: secular: equation. 
is partially diagonal. The energy flow calculations: are- 

‘ a es : - alee 2 ming ae 
simplified since Cr CSE» Czg and C5 ane: identically Zero.» 


Mme remaining constants are 


_ 2 Z 

Cio = (l-n C5 +o C3 

Mme i (i-n-}(c_. + G.. —lc,,) + (i-2n- (iene: 
13 ae 515) yh | ey 

(329) 

_ 2,1/2 

C45 = -n(l-n" ) (Cr _ Cro) 

C25 = Sn (lane Cir - n(Isn*)7/* (1-2 Jc; 4. 

22 vas GENTS) Dips 2a Oe 
#n(i-nP"*c,, — @a(Ton ) “(1-2 C4 

-_ . aaae/2 | 

Bee 7 Con yy - CGg? 


Me resultant velocity and energy flow calculations: are: 
maentical to those for hexagonal symmetry.. The: results:are- 


shown on the next graph for (100) plane propagation: in: Tin. 


mgeeregquired data for calculation. fs giver. below. ° 
Cii= ele7 x idan dynes/em~ | 
Cio = 5-78 i" 
Ci3 = 3.42 u 
C23 =10.31 
Chy = 2.69 
Oe ae oe i 
p = [50 em/em> 
Orpid. 


Be 





The response is gimilar tc thaw [orm Zagreeee 
accidental mode oeceurs.. In adeitaon tre sia meee ee 
angles again indicate that the material is: essentially 
isotropic. For the tetragonal results the: deviation. angles. 


are not reversed in sign... 


Ss 
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VIII... CONCLUSION 


Rereement with published cava: indicates: (nat shee memes 
Preeencved in this paper is 4 valid one and: is= parbicuter 
Weerul in the determination. af energy flow components. Thus 
eG 32S possible to calculate the direction. of! energy fiow for 
Waves propagating in an arbitrary direction... The:results~ 
eso show that, perhaps surprisingly,. there: may be-directions 
m™ preopagation for which all the modes are pure but~ the:-trans— 
mersewwave energy flow vector” and the propagation: vector are 
not colinear. This method can. be extended: to: a. crystal _of 
general symmetry. The effects af crystal stability and multi- 


ple mode response are areas for further investigation. 
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